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ABSTRACT 


A new type of interferometer for use in Fourier spectroscopy has been 
devised at the Aime Cotton Laboratory of the National Center for Scientific 
Research (CNRS), Orsay, France; with this interferometer and with compu- 
tational techniques developed by the Regional Interdisciplinary Center for 
Electronic Calculation (CIRCE) at CNRS, inte rfe rograms comprising as 
many as 10^ samples can now be transformed. The techniques are described, 
and examples of spectra of thorium and holmium, derived from 10^-point 
inte rferograms, are presented. 

(This Technical Memorandum is a translation of an article, 
n Spectroscopie de Fourier avec transformation d'un million de points, n by 
J. Connes et al. , Nouvelle Revue d ! Optique Appliquee, Vol. 1, pp. 3-22, 

1970, published by Masson &: Cie, Paris, France. The abstract above was 
prepared by the translator.) 
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I. INTRODUCTION 


The use of Fourier spectroscopy for problems requiring high re solution 
is now well established and the advantages have been demonstrated. The 
development of the technology is pursued at the Aime Cotton Laboratory and 
a new type of interferometer has been devised; as a consequence of computa- 
tional techniques developed by the Regional Interdisciplinary Center for 
Electronic Calculation (CIRCE) at CNRS, inte rfe rograms comprising as 
many as 10^ samples can now be transformed. 

First, we shall briefly discuss the earlier developments. The first 
generation interferometer (1964-1966) constructed in collaboration with the 
Jet Propulsion Laboratory (Refs . 1-3) provided a re solution of up to 0 . 1 cm 
as a consequence of a carriage movement of 5 cm. It used several novel 
technique s: 

(1) The use of cat's-eye re troreflector s (Ref. 4) in order to elimi- 
nate alignment errors caused by the carriage motion. 

(2) Step-by-step recording using a servo-controlled carriage, the 
signal integration being performed when the carriage was sta- 
tionary. This procedure allowed the path difference to be known 
to a precision approaching 0. 1 nm and used the recording time 
in an optimum fashion by reducing the number of samples N to 
the minimum permitted by the sampling theorem (Ref. 1); N 
therefore becomes equal to the number of spectral elements 

M = Acr/5cr where Act is the spectral interval being studied and 
8(7 is the resolution (or width of the instrumental profile). By 
contrast, workers using continuous scanning must normally use 
N»M. 
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(3) Internal modulation of the interferometer (Ref. 3) (that is, 

modulation of the path-difference itself) at a frequency appropri- 
ate to the detector being used. This technique, which gives, 
approximately, the derivative dl/dA of the interferogram, 
reduces the effect of intensity variations of the source, particu- 
larly those due to atmospheric turbulence. This gives the 
experimenter a free choice of exposure time.* 

The maximum path difference A (= 10 cm) was fixed by the width of 

198 max 

the green line of Hg used for the servo reference, and the maximum step- 
ping rate (5 steps/ second) was determined by the response time of the servo 
loop. The interferograms were recorded on punched paper tape. 

The reductions were performed on an IBM 7040 at the Meudon Obser- 
vatory; the program used allowed 64, 000 points to be transformed (in 
105 min). Despite these limitations, the system was well adapted for the 
study of weak sources such as stars and planets at moderate resolution 
(0 . 1 — 1 cm 1). The planetary spectra we re assembled and published as an 
atla s (Ref . 9) . 


The MK I interferometer showed that the resolving power and the 
instrumental function were rigorously identical to the theoretical values. In 
the meantime, single -mode, frequency-stabilized gas lasers having been 
developed, it became clear that their use as reference sources would allow 
Fourier spectroscopy to be extended to very much higher resolving powers. 


The MK II interferometer, constructed and tested by J. Pinard (1965 
1968) (Refs. 10 and 11) demonstrated the practicability of this extension. 

The carriage displacement was equal to 80 cm, giving A = 160 cm; the 

resolution was 6cr^ 1 /a 
we re : 


max 


6 X 10 ^ cm ^ 


max 

The new techniques developed 


(1) The use of the 632.8 nm line from a single-mode He-Ne laser 
stabilized by the Lamb Dip for a reference source. 


*By contrast, the technique called "rapid scanning" (Ref. 5) connects the 
modulation frequency and the exposure time. Its use on weak sources such 
as stars requires the summation of a great many successive interferograms 
in a digital memory. Although well adapted for low-resolution work, it is 
difficult to use for the kind of problems that interest us. 
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(Z) The use of pressurized oil-bearings for the moving carriage, 

driven by a two-phase linear motor for the servo; the two 

approaches are easily extensible to greater movements, and the 

value of A was chosen only because it seemed adequate for 
max 3 

all conceivable studies. 

The MK II interferometer was conceived for the study of relatively 

narrow spectral regions at high resolution. The recording and reduction 

techniques were essentially the same as for the MK I; since M = N, it follows 

that the maximum number of spectral elements available was M = 64, 000. 

max 

If 6cr = 6 X 10 cm , the spectral interval is reduced to 390 cm~l. This 
instrument, for example, was capable of studying a single vib ration- rotation 
band in the near infrared. 

The results obtained established, on the one hand, that the design reso- 
lution was attained (it is several times better than that of the best grating 
spectrometers) and, on the other hand, that the wave-number accuracy was 
equally improved. In essence, since there is neither scanning of the spec- 
trum nor interpolation between secondary standards, the wave-number scale 

is determined by purely digital operations based upon a single standard line. 

- 1 -4-1 

In a N^O band near 4, 700 cm , an RMS error of 1.5 X 10 cm was 

7 

obtained by comparison to the calculated positions (an accuracy of 3 X 10 ). 

If even greater resolving powers do not seem to be required, then an 
increase in the width of the spectral interval to be studied would be feasible 
and would allow new problems to be attacked at maximum resolution. 

It is clear that Fourier spectroscopy will not have shown its full power 
and the multiplex principle will not be fully exploited until an interferometer 
can, in a single inte rfe rogram, record all the information contained in the 
spectrum of the source to be studied , within the response range of a single 
detector . 

For example, a recording of an emission spectrum ought to permit the 
simultaneous measurement of hyperfine structure (and, eventually, the 
Zeeman structure), isotope shifts, the widths and shapes of the lines, and 
their absolute positions. For sources of limited lifetime, such as radioac- 
tive isotopes (in which we have a great interest), all these data could be 
obtained in a single run. This is impossible by classical techniques. 
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In the study of an astronomical source, one tries to obtain the highest 
resolving power allowed by the available energy in a single observing ses- 
sion; however, the MK I interferometer could only reach 6cr = 0 . 1 cm ^ at 
best and this only within one of the three atmospheric windows accessible 
with a PbS cell (1 to 2.5 pm); hence, three separate runs were needed to cover 
all this region. As for the MK II interferometer, it was never used astro- 
nomically because of the slow recording rate and the narrow spectral inter- 
val at high resolution. 

Such a goal is much more difficult to attain than a greater increase in 
resolving power would be since, as we have shown, such an increase would 
pose no new problems. The problems encountered are of three types; 

(1) The difficulty- -and duration--of the reduction increases with the 
number of points to be transformed and is independent of the 
resolving power. 

(2) The rate of controlled movement of the interferometer must be 
much greater in order to record a given number of data points 
in a given time . However, the accuracy of the system must be 
the same; an error of a single step during the scan results in a 
severe distortion of the spectrum. 

(3) The wider the spectral interval and the greater the resolving 
power, the more the discrimination of the spectral elements 
(which should, in theory, be perfect in a Fourier transforma- 
tion) becomes a practical difficulty. 

It is useful to introduce a quality factor Q for spectra obtained by 
Fourier transformation of an interferogram. This is simply the ratio of the 
total energy 


CO 

B(cr) • do- 

contained in a spectrum to that in the minimum detectable line (absorption 
or emission). With the assumption that the true line width is less than or 
equal to the instrumental width, the strength of such a line will be essentially 
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equal to the RMS value of the random noise in the final spectrum from 

all causes. The quality factor can therefore be written as Q = MS /B rr 

7 mean eff 

where S is the mean intensity in the spectrum, 

mean 7 r 

The only truly fundamental noise of the system is that of the detector. 
However, a number of other sources can exist (variations in source bright- 
ness, amplifier gain or zero level, instabilities in integrators, digitization 
noise, etc.) whose relative importance increases as Q becomes larger. 

If the factor M is a measure of the superiority of a multiplex technique, 
then the factor Q is a measure of the difficulties to be overcome in obtaining 
this gain. The problem which we have attacked here (widening of the spec- 
tral interval at high resolution) required just such an increase in Q for the 
system. The price paid was a greater need for care and a more complex 
system. 

The MK III interferometer, which will be described below, resulted 
from these various requirements. The step-by-step motion and internal 
modulation were retained for their proven advantages, while leaving a much 
wider scope than we had in the earlier systems for a choice of step-length 
and modulation amplitude and being totally secure against external vibration. 

The maximum stepping rate is 50 per second, and 10^ samples can be 
recorded, although this last figure is by no means the limit feasible. At the 
same time, the computational techniques have been developed to transform 
such an inte rfe rogram in 22 min on an IBM 360/75. 
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II. INTERFEROMETER AND SERVO SYSTEM 


A. Optics 

The general layout of the MK III interferometer is similar to that of 
the MK II (Ref. 10, Fig. 2). The maximum carriage movement is 1 meter, 
whence 


A = 2m and 6a = 5X10 ^ cm ^ 
max 

The diameter of the cat's-eyes has been increased from 80 to 160 mm and 
the useful beam diameter from 20 to 80 mm. The single beamsplitter has 
been replaced by two separate plates labelled S (splitter) and M (mixer), 
the thicknesses being identical to within 0. 1 (Jim (Fig. 1). A common sup- 
port and fine flexure adjustments assure their parallelism. The two semi- 
transparent coatings were made simultaneously on both plates, assuring 
their identity and making the system perfectly achromatic. The plates 
themselves are of Infrasil and the coatings are of silicon, the optical thick- 
ness being one -quarte r -wave at 1 . 6 pm; the 4RT product is greater than 
0.75 from 0.8 to 3.5 pm. Other plates are planned for other spectral 
regions . 

The reference line used is the xenon superradiant emission line 
\ = 3.50 pm (Section III- A) . The beam first passes through a linear 

polarizer P (2-mm plate of dichroic calcite [calcium carbonate] cut 
parallel to the optic axis (Ref. 13)). In the two arms of the inter- 
ferometer are two quarter-wave plates and (quartz cut parallel 
to the axis, the thickness being e =136 pm) orientated in such a way 
that the beams are right (D) and left (G) circularly polarized; after 
recombination at M these interfere to give a linear polarization R of 
constant intensity (independent of A) . By contrast, the azimuth n of 
the polarization is a linear function of A and rotates through an angle tt 
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for a change in path difference of X- r ; if is a particular value of A and 
is the corresponding azimuth, one finds that 

tt( A - A q ) 


Fixing the interferometer at a particular A is equivalent to requiring that 
ex take a given value. 

The measurement of a is reduced to the measurement of the phase of 

a modulated signal . A half-wave plate D, which rotates at constant speed 

(N turns/second), rotates the transmitted plane of polarization at 2N turns/ 

second. After passing through a fixed linear polarizer A (identical to P) 

s 

the flux is sinusoidally modulated in intensity at a frequency of 4N Hz*. A 

detector H (an uncooled InAs photovoltaic cell) gives a signal F at a fre- 
s s 

quency of 4N Hz which acts as a phase-modulated carrier wave for the 
path-difference variations in A . If the zero of phase occurs when A = A^, 
then we may write 


F 

s 


= I sin 2 7 r 


4Nt + 


A ~ A o\ 

X r / 


B . Error Signal 

From F g one seeks to produce an error signal , that is, a dc volt- 
age E proportional to the difference (A - A^); furthermore, it is desir- 
able to reduce the residual ripple of the 4N Hz carrier wave as much as 
pos sible . 


*This is the same as saying that the half-wave plate increases the fre- 
quency of one of the circular polarizations by 2N Hz and reduces the fre- 
quency of the other by the same amount; the detector H s detects the beats 
between these two frequencies. 
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The signal F is first of all shaped (Figs. 2 and 3) and results in a 

square wave of the same phase and of unit amplitude C g . A synchronous 

demodulator multiplies C with a square-wave reference R having the 

same frequency and an adjustable phase whose production is discussed 

later. The product P = C • R can be considered as a train of positive 

or negative pulses of unit amplitude, of which the width depends on the 

relative phase of R and C . The mean value is zero if C and R are in 
r s s s s 

quadrature, equals + 1 if they are in phase, and equals - 1 if they are in 

opposition. This mean value can be used as an error signal; its variation 

as a function of A (for a given phase of R g ) is shown in Fig. 4: it is linear 

over a range of ±\^/4 variation in A around the value Aq at which C g 

and R are in quadrature, the nominal operating point, 
s 

However, at this condition there is an unwanted residual ripple 

superposed which has its greatest amplitude just at the position that A = A^; 

this ripple is a square wave of frequency 8N and of unit amplitude. In 

order to suppress it, a second detector is used which receives half of 

the reference beam by way of a semir eflecting mirror L and a second 

analyzer A , which is set at 45 deg with respect to A (Fig. 1). A second 
c s 

signal F , in quadrature with F , is thereby generated; after shaping, it 
gives C^, in quadrature with C . A second synchronous demodulator gen- 
erates the product P = C • R , where R is a second reference signal 
^ c c c c 6 

automatically in quadrature with R . The sum E = P + P is taken and 

s c s 

used as the error signal. Figure 2 shows that E is composed of a train of 
rectangular pulses, identical and equidistant, of constant height, of fre- 
quency 16N, of positive or negative polarity depending on the sign of 
( A - Aq) and of width 


T =|n.| a - A 0 |A r 

When A = Aq, E = 0 and the carrier is totally suppressed. For 
A ^ Aq, the mean value of the pulse s - -which is the required information- - 
is once more given by Fig. 4; a residual of the carrier remains but at a 
frequency 16N; the passband of the servo system is, however, funda- 
mentally much less than 16N. It is equivalent to saying that if, at a 
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certain instant, a step-function error A is imposed upon the rest value A^, 
the information necessary to make the correction will be unknown until 
after an error pulse has been received; that is, after a delay which is at 

most 1/16N. With the value of N = 1000 turns/second finally adopted, this 

-4 

delay equals 0.62 X 10 seconds; it is nearly negligible compared to 
the mechanical response time of the system. 

C. Damping Signal 

In order to obtain an optimum servo feedback, it is useful to have a 

damping signal ; that is, a second dc voltage E 1 proportional to dA/dt. 

Noticing that if the path difference changes at constant speed, following a 

law ( A - A n ) = Vt, the phase of the signal F will vary linearly with time; 
u s 

that is, the frequency of the carrier is changed; it becomes equal to 
4N^ = 4N + V/\ (single side-band modulation). 

Figures 5 and 6 explain how E* is produced. The signals C g and 

come from the infrared detectors at a frequency 4N^ (variable with V) and 

are used to generate two pulse trains I and 1^ at twice the frequency 8N^ 

(by differentiation and full-wave rectification) displaced one from the other 

by a half-period. Each pulse-train unlatches a single- shot in which the 

excited state has a fixed delay 6=1/ 16N, which is exactly equal to half the 

interval between successive pulses 1/16N when V = 0 and N = N. 

v v 

If N < N, 0 > 1/16N , and if N > N, 0 < 1/16N . The output volt- 

V V V V 

ages of the two single-shots are and M g ; the quiescent state corre- 
sponds to a voltage of - 1 and the excited state to +1. The summation 

results in a train of rectangular, equidistant pulses of constant 
amplitude and of frequency 16N • The pulses are positive if N^< N and 
negative if > N; their width is 


T* = I 1/16N - 1/16N I 

1 v 1 

If IN - N|«N, T' « 1/64 • v/n 2 x. . 

' v ' r 
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The mean value of the train is therefore proportional to V, and when 
V = 0, T ! = 0. This mean value constitutes the desired damping signal E ! ; 
the carrier ripple, of frequency 16N, is again zero at equilibrium (V = 0) o 

Up till now, we have assumed that the rotational frequency N was 
perfectly constant. It is easy to compensate for the effect of its variation 
if the reference signals and R g are generated from the same device. 

The frequency is 4N (independent of V) and the pulse train thus produced* 
is subtracted from the pulse train discussed above (by sending it to the 
negative input of a differential amplifier). The servo system has two 
parallel loops. The low-frequency components of the error signal (roughly 
0 to 50 Hz) are sent to a linear motor (see Section III-B) which controls the 
position of the moving carriage (mass = 1Z kg). The high frequencies 
(above Z0 Hz) are sent to a pair of piezoelectric elements which carry the 
small secondary mirrors of the cat's-eyes (mass = 1 g) for which the max- 
imum displacement is ±20 pm; a step function in path difference A of ±80 pm 
can therefore be made without moving the carriage. 

D. Reference Signals R g and R^ 

An air-driven turbine with a hollow axis carries both the half-wave 
plate D and a glass disk upon which are 400 opaque radial sectors. Two 
small incandescent lamps and mounted diametrically opposite each 
other illuminate two silicon detectors S^ and S^ (Fig. 7); a second, fixed, 
disk Vg is mounted coaxially to V^. The sum of the two signals from S^ 
and S^ is used to eliminate eccentricity errors. 

The signal K thus produced has a frequency of 400N, being 400 kHz if 

N = 1000 turns/ second. The reference signals R g and R^ are generated by 

dividing K by a factor of 100. The usefulness of this approach lies in the 

ability to alter the phase of R and R by digital methods; the phase can 

s c 

take any one of 100 values incrementable in steps of exactly 1/100 of a 


*The mean value of this pulse train varies with N; it can be used to servo- 
control the rotational frequency of the half-wave plate. To date, a simple 
control of the compressed air supply giving a frequency stability of about 
1 % has sufficed. 
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period. Consequently the operating point A = A n (cor re sponding to C and 

u s 

C being in quadrature with R and R respectively) can be shifted by 
c sc 

hundredths of a reference fringe. The smallest sample step is therefore 
equal to X^/100 (35 nm). 

Figures 7 and 8 explain the technique. The signal K is squared and 
differentiated. Two 400-kHz pulse trains are thereon switched into the * 
lines A and B by inverters; the pulses A and B are displaced by a half- 
period and cannot coincide. The pulses A are sent directly to the positive 
(counting) input of an up-down counter, providing a division by 100; the 
output pulses actuate two switches which produce the square waves R and 
R c> automatically in quadrature. The pulses B can be sent by the gates 
P + and P to either the positive or the negative inputs of the same counter. 

If P + and P are closed, the frequency of R and R is exactly 1/100 
— sc 

of that of K (400N/100 = 4N) ; their phases are constant and the inter- 
ferometer carriage does not move. If P + passes a single B pulse, the 

phases of R_ and R_ advance by 1/ 100 of a period and A increases by an 
s c 

amount equal to X / 100. If P is used, an equal and opposite decrease 
in A is obtained. If P + is held open, the frequencies of R g and R^ become 
8N; if P , they fall to zero; the interferometer carriage thereupon goes 
into a constant velocity condition, forward or backward. The correspond- 
ing velocity in A is ±4NX^ ( 14 mm / s for N = 1 kHz and X = 3.5 pm) ; this 
is the maximum (slew) speed that is permitted by the system. 

E. Path Difference Scan Program 

By means of an appropriate programming method for P_^_ and P , any 
path difference scan program Aft) which satisfies the foregoing restric- 
tions (minimum step 35.0 nm, maximum velocity 14 mm/s) can be obtained. 
In particular, the internal modulation itself can be included into the pro- 
gram and the amplitude controlled by the servo system, which was not the 
case with the earlier systems. 

The program A (r) actually used is represented by Fig. 9. The oper- 
ating cycle is made up of two parts. During the interval AB the carriage 
is at rest and the internal modulation is performed by the piezoelectric 
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elements. Their peak-to-peak amplitude is q\ /100 (where q is integral and 

3 r 

< 10 )*; the frequency N ! is a submultiple of N (we generally use N ! = 250 Hz). 

The resting time AB (and the integration time for the signal coming from the 
source) is an integral multiple of the internal modulation period 1/N 1 . Very 
short integration times can therefore be used without error (Ref. 2, p. 126). 

During the interval BC, A increases by P\ / 1 00 (where p is integral and 

4 r 

<10 ), the basic step-length, ** and the integrator is reset to zero. The step 

is produced by a movement of the carriage; however, because of the inertia, it 
follows the scan program A (t) only after a delay of a few milliseconds, and the 
piezoelectric elements are used to compensate for this automatically. The 
maximum scan rate is 200 steps/ second. 

For problems of a type somewhat different from those envisaged here 
(intense sources of short lifetime such as plasmas, for example) one might 
try for a much higher stepping-rate (the number of samples N would not neces- 
sarily be greater). The only fundamental limitation would seem to be the 
response time of the infrared detector being used. However, under these 
conditions, there would no longer be any reason to keep the step-by-step scan 
and the internal modulation; the technique of "rapid scanning" (Ref. 5) would 
seem to be more appropriate. Moreover, there is no discontinuity between 
the approaches; when the integration time AB becomes equal to a single per- 
iod 1/N 1 , the recording rate with internal modulation is at its maximum; the 
(much greater) recording rate becomes, in effect, a rapid scan. Our servo 
arrangement can work equally well in such a mode but we should need a much 


*It should be remembered (Refs. 2 and 3) that the amplitude of the internal 
modulation must equal \ ean^' w ^ ere ^mean ” mean andcr mean the 
central wave number of ffi| a spectral interval under study, and that a square - 
wave modulation is optimum. However, the advantage over a sinusoidal 
modulation is small and is not enough to make the effort needed to generate a 
good square worthwhile. In the near infrared (where we plan to work) the dif- 
ference is negligible. For instance, if \ mean = 10 fim, the modulation ampli- 
tude is 2 .5 pm and the rise time 0 1 = 3,5 X 10‘^ s is less than a tenth of a 
modulation period (for N = 1 kHz and N' = 250 Hz). 

**The rules for determining the step-length have been given elsewhere 
(Ref. 12). It may be recalled that it is determined by the limits c r j and cr2 of 
the spectral region under study; to a first approximation it is of the order of 
1 / 2 ( o*2 - 0 * 1 ) = 1/2 Ac r . 
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faster incremental recording system than the present one, which is limited 
to 50 measurements/ second (see Section III-E). 

F . Advantages of the System 

The first advantage is the smallness of the minimum step X^/100 = 

35.0 nm, which allows a much more flexible choice of sampling interval than 
the earlier systems; one may therefore use only the theoretically minimum 
number of samples N as allowed by the sampling theorem (which is essential 
when the number is very large). Similarly, the controllable amplitude of the 
internal modulation allows optimum usage for any spectral region.* 

A second, less obvious, advantage is the following: in spite of the 

smallness of the basic step-length, the functioning of the system is much less 
susceptible to vibration than the MK I or II because we do not count interfer- 
ence fringe s . The counting depends only on the 400 -kHz signal, for which the 
amplitude and frequency are stable to about 1 % and the signal -to -noise ratio 
is very high. Intensity fluctuations in the superradiant tube can in no way 
introduce a counting error. 

The response range** of the servo is rigorously linear over a change in 
A of ± \^/4 (±0,88 pm). For the MK I and II interferometers, the error sig- 
nal was generated starting from the intensity of sinusoidal fringes; with the 
632.8-nm line from a He-Ne laser, the useful response range (roughly 
linear) was about ±0.1 pm. 

These improvements are directly responsible for the improvement in 
the servo system bandwidth and the gain in speed and security of the entire 
system; they have allowed us to go from 5 points/ second to 200 points/ 
second and from an N of 60, 000 to N = 10^. 


*The MK II interferometer necessarily used an amplitude that was a multi- 
ple of the reference wavelength. 

**This range can be increased without limit should the need be felt. Starting 
from the already available signals C s , C c , R c , and R s , a second error sig- 
nal Ei, in quadrature with E, can easily be generated (Fig. 4). The devel- 
opment of E and El starting from a position A q defines the value of A without 
ambiguity; a reversible period counter on E and Ei would permit an error 
signal to be produced that cancels only once (at A =Aq). 
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Fig. 1. Simplified schematic of the optical system. Only the reference 

beam (10 mm diameter) is shown; it travels down the center of the 
main source beam (80 mm diameter). The interferometer, on the 
one hand, and the super radiant tube, the turbine and the detectors 
on the other hand, are mounted on separate baseplates. 

NOTE ; OEIL DE CHAT = cat’s eye 
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Fig. 2. Method of generation of the error signal E = P g + P c> starting 
from the phase-modulated signals F g and F c and the reference 
signals R g and R c . The potentiometer P allows a dc voltage to be 
added to the error signal and a small change in A, equal to a frac- 
tion of the elementary step-length X.^/100, to be produced; it may 
be directly calibrated in terms of path-difference. It is used to 
place the first sample in the interferogr am exactly at the point 
A = 0. 

NOTE : ECRETEUR = shaper 

CERAMIQUE = piezo -electric element 
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1/4 N 



Fig. 3. Waveforms of the signals from which the error signal E = P g + P c 

is generated. The figure illustrates the condition that A - Aq = 

\ / 1 6 ; that is, C c and C are out of phase by 1 / 1 6 of the period 
r sc 

with respect to the quadrature with R s and R c (which would give 
E = 0). 
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Fig. 4. Variation of the error signal E as a function of A around the oper- 
ating point Aq. The phase of the reference signals R g and R c is 
presumed constant. The function is periodic with a period it is 
linear within a range of X- r /2. The error signal has the correct 
sense (and the restoring force is towards Aq) over a range of \ r . 

A second error signal in quadrature with E is represented by Ej 
(see footnote, page 13). 



Fig. 5. Generation of the servo feedback signal 

NOTE ; DOUBLER = frequency doubler 

MONOVIBRATEUR = single- shot multivibrator 
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Fig. 6. Waveforms of the signals producing the damping signal E* = 
M c + M g . The figure shows the case wherein N v < N. 



Fig. 7. 


Generation of the reference signals R g 


NOTE : LAME DEMI-ONDE 

DISQUE FIXE 
DISQUE TOURNANT 
ECRETEUR 
COMMANDE PORTE 
OU 
ET 

COMPTEUR REVERSIBLE 


and R c at a frequency of 4N. 
= half-wave plate 
= fixed disk 
= rotating disk 
= shaper 

= gate-open pulse 
= (logical) OR 
= (logical) AND 
= up- down counter 
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r s < {2) 2ZZ 1 

J3) 

Fig. 8. Reference signals R g and R c . Starting from the signal K, at a fre- 
quency of 400 N, two pulse-trains A and B are generated. The 
leading and trailing edges of the square-wave R g (1) leaving the up- 
down counter are defined by the pulses A. If the gate P allows an 
extra pulse to pass, the phase of R g is retarded by l/lOO of the 
period (R (2)); if it is the gate P , the phase is advanced by the 
same amount (R g (3)). The waveform R c (not shown) stays in 
quadrature with R g . The figure is not to scale, the ratio of 100 
between K and R g not being maintained. 
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Alt) 



Fig. 9. The function A(t). The finite (error) width of the elementary step is 
neglected in this figure. The symbol A ! (t) represents the actual 
carriage motion; the displacements of the small mirrors mounted on 
the piezo-electric elements are given by A - A*. 

NOTE ; DEPLACEMENT DU CHARIOT = carriage motion 
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III. BRIEF DESCRIPTION OF THE SYSTEM ELEMENTS 


Only the novel items are discussed here; a more detailed description 
will be given in subsequent publications. 

A. Super radiant Xenon Source 

This is a tube 1 m long and 3 mm in diameter (just enough to pass the 
TEMqq axial mode at 3.5 pm); the useful length is doubled by virtue of a 
fully reflecting mirror at one end. Filled either with xenon or a Xe-He 
mixture and dc excited, this tube emits an intense 3.5-pm line along its 
axis . 

-3 -1 //t 

The calculated Doppler width is of the order of 4 x 10 cm (for 
400 K); the effective line width is much less than this because of the super- 
radiance. In fact, neither the width nor the profile have been measured 
since these parameters do not intervene in the use of the line for a refer- 
ence; a small reduction in depth of modulation is noticed at A = 2 m; this 
reduction has no effect on the phase modulation system. 

The position of the line is not defined by a cavity* and therefore con- 
stitutes a natural standard. The principal source of instability appears to 
be variations in pressure; the pressure is maintained at the saturation 
vapor pressure of xenon at a well-determined temperature (about 77 K) 

using a reservoir of xenon submerged in liquid nitrogen boiling under a 

2 

pressure of about 2 kg/ cm . To date, only natural xenon has been used and 

1 36 

the absolute wavelength \ has not been measured. Eventually, xenon 
will be used. 

This tube is a convenient reference source for Fourier spectroscopy, 
much simpler than a single-mode servo-controlled laser. Because of the 
greater wavelength, the servo control is more certain; furthermore, mate- 
rials which are not transparent in the visible (germanium or silicon) or which 
cannot be polished to a sufficient precision (sodium chloride) become usable 
as beamsplitters. 


*It is desirable to suppress light reflected back to the super radiant tube by 
the various optical elements which can result in laser-type action. In our 
equipment, the quarter-wave plate Q placed in front of the polarizer P 
(Fig. 1) forms an optical isolator. 
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Linear Motor 


B . 

The linear servomotor for the MK I interferometer was a simple 
solenoid in the radial field of an annular magnet; the range (5 cm) was lim- 
ited by the length of the solenoid. That of the MK II interferometer (Refs. 10 
and 11) was a two-phase asynchronous type and the maximum available thrust 
was only a fraction of a newton force. This type of motor has the inconveni- 
ence (important for a system to be operated in a vacuum) of a significant 
power dissipation (even when the thrust demand is zero) because of the ac 
polarizing current. 

A dc linear motor, with a noncontacting commutator for detection of 
the field polarity by means of Hall effect detectors, has been employed for 
the MK III Interferometer (Fig. 10). The maximum thrust available is 
40 newtons force for a dissipation of 100 W (this value being the limit of the 
dc power amplifiers). In fact, the carriage acceleration is found to be lim- 
ited by mechanical resonances and the maximum thrust is never used. We 
now think that a linear dc motor of a more classical type (long pole, without 
commutation) could suffice. 

C . Piezoelectric Elements 

These are the high-frequency elements of the servo system. Fig- 
ure 11 shows the layout. They drive the small secondary mirrors of the cat's- 

eyes. The elements are enclosed in copper tubes; a silicone damping oil 

Z 

(viscosity 0.1 m / s) provides simultaneous mechanical damping and thermal 
dissipation. The whole arrangement can operate in vacuo . The drive for 
the electrodes comes from high-voltage transformers. The maximum allow- 
able voltage between the internal and external electrodes is 2000 V and pro- 
duces the 20-fJim displacement cited earlier. 

D . Rotation of the Half-Wave Plate 

A commercial grinder bearing with a hollow axis and driven by com- 
pressed air is used for the rotation of the half-wave plate. The noise level 
and the vibration are tolerable and cause no particular problem. However, 
the lubrication (and lifetime) of the ball-races is unsatisfactory: the use of a 
rotary air-bearing should resolve this difficulty. 
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E . Measurement and Recording Arrangement 

Figure 12 shows a simplified schematic of the system arrangement. 
The signal to be measured is demodulated by a 250 -Hz square wave, 
of adjustable phase, synchronous with the internal modulator. Integra- 
tion is performed by a voltage-to-f requency converter followed by a 

5 decimal digit up-down counter; the samples are stored and a new 
integration begun during the writing-time on the magnetic tape of the 
previous sample. Each measurement involves the recording of 6 binary 
coded decimal (BCD) characters, one of which is a flag character indi- 
cating the gain of a manually operated 4-position attenuator (gain of 1, 

10, 100, 1000); its purpose is to augment the dynamic range of the sys- 

tem (Ref. 2). 

The total intensity of the source under study is determined by a 
detector R 1 , and a voltage proportional to it is sent to a conventional 
potentiometric recorder P, which gives a permanent indication of it. A 
second potentiometer, coupled to the first, is introduced into the feed- 
back loop of a variable -gain amplifier; source intensity fluctuations are 
thereby compensated (on the assumption that they are achromatic, of 
course) . 

The recorder is of the incremental type* (400 steps/ second). When 

6 characters are recorded, taking dead-time into account, the maximum 
rate is 50 points/ second . 

F. Adjustments 


The system of phase modulation is, by its very nature, insensitive to 
intensity variations in the reference signal; it requires neither a delicate 
balance between two different detectors, nor manipulation of precise 


*Step-by-step recording allows stops of any length to be made during a 
recording. Thus magnetic tapes directly readable by a computer can be 
made, which must necessarily be composed of data records separated by 
gaps. If the interferometer ran continuously, a buffer store would be 
necessary. 
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waveforms by very linear circuits;* the analog part is reduced to an absolute 
minimum (preamplifiers) and all subsequent operations are digital. How- 
ever, it must be expected that any system of path difference measurement 
not operating by simple counting of whole fringes (and thereby allowing an 
estimate of the fractional error to be made) must introduce a periodic error 
of period ours is no exception. 

Imperfections of the optical system have two different consequences: 

(1) The calcite dichroic linear polarizers are of very high quality 
(Ref. 13), but by contrast the thickness of the quarter-wave 
plates is accurate only to about 1%. To try to get better preci- 
sion would be pointless, because the number of oblique reflec- 
tions in both the interior and the exterior of the interferometer, 
not to mention the oblique transmission through the plates S and 
M, causes the interfering vibrations to be deformed so that the 
polarization becomes elliptical instead of circular. 

It is fortunately possible to compensate for all these errors 
without measuring them individually, by making the incidence 
on the quarter-wave plates and slightly oblique. This 

compensation is very delicate; it requires, in effect, a simul- 
taneous 3-parameter adjustment (azimuth of the quarter-wave 
plate in its plane, the angle of incidence and the plane of inci- 
dence). It is performed by occulting one of the beams in the 
interferometer and orientating the quarter -wave plate in the other 
beam in such a way as to annul the modulated signal received by 

either of the detectors H or H ; the incident vibration on the 

s c 

rotating half-wave plate is then necessarily circular. 


*These qualities, which are essential for a precise measure of A , are a 
property of the mechanical system which produces the modulation (rotating 
halfwave plate D mounted rigidly to the disk Va) * Of course, it is possible 
to perform the equivalent operation with electro-optic crystals driven by 
suitable sinusoidal voltages. But in order to obtain a rigorous relationship 
between the phases of the reference signals of frequency 4N and 400N, an 
extreme precision in the waveforms and amplitudes of the applied signals, 
not to mention the driving conditions and the crystal temperature control, 
would be necessary. An alternative measurement system, somewhat like 
ours (Ref. 14), in which the carrier is obtained by beats between the Zeeman 
components of a laser line, has the advantage of not requiring a rotating ele- 
ment. But it no longer allows the fine, accurate, and rapid interpolation 
between whole fringes which is the essential property of our arrangement. 


24 


JPL Technical Memorandum 33-525 



This empirical adjustment is convenient only if it is possible to 
follow the phase and the amplitude of the modulated signal; there- 
fore, the signal (before squaring) is sent to two synchronous 
demodulators fed by the reference signals R g and R^. Their dc 
outputs are fed to the X and Y inputs of an oscilloscope and one 
attempts to center the spot. The accuracy of the adjustment is 

limited by the signal -to -noi se ratio of the reference signal; an 

-4 

ellipticity of less than 1*5 x 10 has been obtained, which cor- 
responds to a periodic error e of the order of 0 . 5 nm in the path 
difference. Such an error must introduce ghosts into a spectrum 
calculated by Fourier transformation* for which the spacing from 
the main line equals cr^ = l/\^ and the intensity is calculated as 

2tt€ /\. The measured ghost intensity for a visible line (k = 

_ 3 

632.8 nmjwas 5 x 10 ; since it is proportional to 1 / X it is much 

reduced in the infrared. This intensity can appear relatively 
large; but (because of the use of a phase-sensitive Fourier 
Transform) the resultant shape of the ghosts is asymetric and 
they cannot be confused with normal lines. 

(2) The adjustment of the incidence angle on the rotating half-wave 
plate D is significantly more delicate. But the error in A intro- 
duced by the error in the optical thickness of D is time -pe riodic 
(of period 1/N). Since the internal modulation period 1/N* and 
the integration time is an integral multiple of 1/N, the average 
residual effect on each inte rfe rogram sample is nil; in other 
words, only the form of the internal modulation program is 
(very slightly) changed and is unimportant. 


*Since the periodic error is stable, one might imagine measuring 
compensating in order to reduce the ghosts even more. One can, 
add a dc voltage, which is a function of the fractional error in the 
fringes, to the error signal. One can also introduce a correction 
computer reduction program; these two approaches have not yet 


it and 

for instance, 
reference 
into the 
been tried . 
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Fig. 10. Direct current linear motor. The interferometer carriage carries 
a coil B and a Hall effect probe H. Two series of stationary 
permanent magnets are mounted on each side of the coil in the 
direction of motion. The magnetic field is periodic; the size of the 
pole-pieces and of the coil are such that the force generated with a 
constant current is nearly sinusoidal. The probe senses the zeros 
of the field and causes the coil current to be reversed; the curve 
showing the force as a function of displacement is therefore the 
rectified sine-wave F^ . Three identical combinations (coil and 
probe) are used, simultaneously, displaced by l/3 of a period, 
giving the forces Fj , F^, and F^. The resultant force F^ is 
thereby only slightly modulated. 
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Support isolant 


Joint dfetanchftite 



Fig. 11. Mounting of the piezo -electric elements 


NOTE : JOINT D*ETAN CHEITE ELASTIQUE 
SUPPORT ISOLANT 
POLARIZATION 
GAINE EN CUIVRE 
HUILE D'AMOR TISSEMENT 
CERAMIQUE 
MIROIR 


watertight, elastic, joint 
insulated support 
power leads 
copper sleeve 
damping oil 
piezo-electric tube 
mirror 
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Fig. 12. Schematic of the measuring and recording system (analog and 
digital parts) 


= control laser 


NOTE ; 

SOURCE DE REFERENCE 

SOURCE ETUDIEE 

MESURE D'INTENSITE 

SYSTEME PAS- A- PAS 

ATTENUATEUR MANUEL 

ATTENUATEUR ASSERVI 

POTENTIOMETRE ENREGISTREUR 

CONVER TISSEUR TENSION- 
FREQUENCE 

R.A.Z. 

INDICATION DE LA GAMME 

COMPTEUR 

CODEUR 

ORDRE D'ENREGISTREMENT 


- source being measured 
= total input power monitor 
= step-and-lock logic 
= manual attenuator 
= servo-controlled attenuator 
= potentiometric strip-chart recorder 

= voltage-frequency converter 
= reset to zero 
= gain indication 
= counter 

= encoder (coupler) 

= write command 
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IV. METHODS OF REDUCTION 


The problems of the redaction process can be classified into two 
categorie s: 

(1) The calculation of a clean Fourier Transform (FT). 

(2) Various operations performed on the recorded inte rfe rogram 
before calculation of the FT; interpolation of the output spectrum 
between the calculated primary points after the FT; plotting the 
spectra and their exploitation. 

Following the programming of the mathematical principles involved in 
the treatment of the inte rfe rogram and of the spectrum, the principal diffi- 
culty lies in the number of data points to be manipulated, which is always 
greater than the size of the core of a computer. We have transformed 
60,000 points with an IBM 7040 having a 32K word* core. We now have a 
360/75 having 256K words but we need to transform 10^ points. 

Now, the fast FT method of Cooley and Tukey is easily applied when 
all the computation can be done in core. However, when the number of 
points exceeds a limit N^ f depending on the computer being used, auxiliary 
stores (tapes or disks) must be used, for which the access time is much 
greater than that for the core. It has been necessary, therefore, to gene- 
rate two programs for transforming a number of points which, in one case, 
is less than N^ and in the other, is greater. 

A. Calculation of the Fourier Transform 

We must calculate the sine FT of an odd function sampled by N points 

between A = 0 and A . Let us put 

max r 

P\ , 

1 _ r _ 1_ 

100 " 2Acr 


*1K = 2 10 = 1024. 
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where h is the distance between sample points (Fig. 13). Starting from 
the N samples of the interferogram I (A), one calculates N samples of the 
spectrum S (cr) separated by 6cr by the series 

k=N 

S (c Tj + p6cr) = 1^. J cos [2tt(cTj + p6cr) (k-1) h] 

k=l 


1, Classical method . The classical algorithm using this formula 
leads, for the calculation of the N points in the spectrum, to the product of 

a square matrix of rank N by a column matrix of N terms, wherein the 

2 

calculation time is proportional to N and the coefficient of proportionality 
is dependent on the computer used and the method of generating the cosines. 
The best times realized follow the formulae 

T = 0. 8 X 10 -3 N 2 for an IBM 704 


and 


T = 0. 3 X 10‘ 3 N 2 for a 7040 

with T being expressed in seconds (Fig. 14 and Table 1). 

In practice, the largest scale on which this algorithm was used was 
for the calculation of spectra of Jupiter (1965) with N = 12, 000 in 1 2 h. 

2. Cooley- Tukey method . In its original form, this method permits 
the calculation of the Fourier Transform of any kind of complex function. 
Thus, in our notation, = 2N complex samples are transformed in order 

to obtain complex values of the spectrum spread out in the interval 
2Ac r = l/h. 

The general idea of the method was expressed at the beginning of the 
fundamental paper by Cooley and Tukey (Ref. 24): "the square matrix 
Nr X N-j-^ can be considered as the product of m partially-filled matrices, 
m being equal to log^ N^. M 

a. Complex functions. In this case, the necessary core size is 2N-j- ) 
words, which implies reserving 4N words for the 2N real samples taken 


30 


JPL Technical Memorandum 33-525 



between -A and +A * and ZN imaginary, null values. Let H be the 

max max 7 

size of the memory of the computer under consideration; H is always a 
power of Z. Since a portion is necessarily taken up by the system monitor 
and the program, one can, at most, use H/ Z words for storing the data; 
therefore with a 7 040 with a 3ZK core one can, at most, handle an inter- 
ferogram of 


N = — = 4K = 4096 words 
-U Z A 4 

With a 360/75 with a Z56K core, N becomes 3ZK. 

L 

On the other hand, the results of the FT are not in order and the 

2 

elementary sorting by inversion of the bits adds a term in to the cal- 
culation time which is quite negligible for small but becomes dominant 

for large N . An improvement of the sorting method has required the use 

u 2 
of an extra words in the core but has suppressed the term in N.^. Two 

equivalent programs, an FFT for 4K words and an FFT for 3ZK words, 

have been written for a 7040 and a 360/75 and have resulted in the following 

calculation times, in seconds, respectively: 

T = 9. 428 X 10" 3 N d log N D + 0. 166 X 10" 2 N D 


and 


T = 0. 152 X 10" 4 N d log ^ N d + 0. 145 X 10" 3 N D 

It should be noted that the terms in log^ and in are com- 
parable and that the change from a 7 040 to a 360/75 gave a speed gain of 
Z3 (see Fig. 14). 

b. Real functions . It is possible to convert the calculation of the FT 
of a real function I 1 , sampled at points, to the calculation of the FT of a 
complex function I u sampled at N^Z points (Ref. Z5). The real part of I" 


*It may be shown that this is equivalent to writing (in the core array 
reserved to the real values) N samples taken between 0 and ^ ma x followed 
by N zeros. 
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comprises the samples of even position and the imaginary part those of odd 
position. Afterwards, it is necessary to re-order either the first N^^/2 real 
samples or the last Nj^/2 of the spectrum of I' , depending on the position of 
the spectrum to be studied in the free spectral interval l/h, starting from 
the complex values of the spectrum of I". By applying this technique, 

we have written a 64K FFT program for the 360/75 which has two advantages 
over the earlier methods: it uses only 2N words of memory for storing the 

data, thereby allowing us to handle = 64K within the central processing 
unit, and the execution time for a given N is divided by 1. 5. The time 
formula (in seconds) is 

T = 0. 158 X 10" 4 N log ^ N + 0. 206 X 10" 3 N D 

c. Extension of the Cooley-Tukey method for the case of N > N^. In 
order to handle the case of N > 4K with the 7040 or the case of N > 64K with 
the 360/75, auxiliary storage must be used. On the other hand, at each 
stage of the Cooley-Tukey method, the numbers to be handled in succession 
are not stored sequentially in the memory, so that it is not possible to use 
the all-encompassing prior programs directly. A different program was 
written for the 7040. It used four magnetic tapes and the final sorting was 
not done by bit-inversion but by a formula from which the position of each 
output sample could be calculated. The running time is given by the 
expression 


T = 0. 228 X 10' 2 N d log ^ N D + 0. 747 X 1 0~ 2 N D 

If N = 64K, T = 105 min. It would be possible to make a significant 
reduction in the coefficient of the first term but this was not done because a 
much more powerful computer became available. 

For the 360/75, we have written a Fourier Transform program for 
1000K points using the so-called ’’decimation in time” technique (Ref. 26), 
which allows an FT of a real function, sampled at = 2^ points, to be 
calculated on a computer having only N^/l2 words of fast-access memory. 

In the first stage, the samples are sorted and separated into 
16 blocks of length Nj_y/l6. These blocks are stored on direct-access 
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locations on a magnetic disk. In the second stage, the Fourier Transform 
of each of these blocks is calculated, in succession, with a 64K FFT pro- 
gram. The third and final part comprises a series of four stages of inter- 
combining the numbers thus produced. This allows the reconstruction of the 
FT of the whole. The spectrum samples (N^y^ for the real part and N^/ 2 
for the imaginary part) are thereby ordered sequentially on two locations of 
the disk. 

The total execution time is the sum of the execution times (Tj, T^> 
and T^) of the three parts: 

Tj = 0. 562 X 10" 4 N d 

T 2 = 0. 151 X 10“ 4 N log 2 N + 0. 206 X 10' 3 N D 

T 3 = 0. 336 X 10' 3 N d 


whence 

T = 0. 151 X 10" 4 N log 2 N + 0. 119 X 10" 3 N 
giving for 

N = 1024 K, T = 22 min, 08 s 

More than half the time is used in the manipulation of the numbers, 
which explains why, starting from N = 32K, the time increases proportion- 
ally to N. 

The size and number of locations has been specifically tailored to the 
case of N = 1024K and is particularly mismatched for N < 64K, but then one 
may use the normal 64K FFT. A factor of 3. 4 is lost in execution time in 
going from a 64K in-core calculation to that using disks, but the core usage 
is 5N/ 2 in the first case and is only 5N/32 in the second. 
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Examination of Fig. 14 shows the progress that has been made, thanks 
to the evolution both of computers and of computing techniques. Fourier 
Transform spectroscopy is no longer classified in the group of techniques 
that are gross consumers of machine time. 

These figures are not optimum. The first improvement that can be 
made is to adapt the FFT algorithm to the particular case of a real function 
which is partly even and partly odd. It is reasonable to expect that this would 
result in less memory being used and the computation time being reduced. 

B . Auxiliary Calculations 

The calculation of the Fourier Transform itself is only one stage in the 
complete treatment which the interfe rogram must undergo between recording 
it on magnetic tape and plotting the final spectrum. As was discussed ear- 
lier, the inte rfe rogram is sampled at the minimum possible number of 
points N. Consequently, the spectrum obtained by the Fourier Transform is 
also sampled at N points, N always being a minimum. This signifies that 
there is only 1 point per spectral element; a linear interpolation gives only a 
rough approximation to the true spectrum and a precise interpolation scheme 
(by convolution) is indispensible (Ref. 1). 

The computation is performed as follows: the 1024K primary spectral 
points S(ct) resulting from the FT are stored on disk. The search for the 
maximum value of S takes 25 s; the normalization of all the spectrum and 
the writing of it onto magnetic tape in IK records takes 86 s; the calculation 
of the requisite interpolation function f((r) and various necessary initializa- 
tions uses another 25 s. At this point, the interpolation proper begins; one 
may calculate either the entire spectrum or an arbitrary portion; the com- 
pie te calculation, employing 5 secondary points per primary point, giving 
5 X 10 24K~ 5 X 10^ points, takes 22 min, 10 s . * Finally, the writing of the 
magnetic tape for the Calcomp 563 plotter is performed in 10 min.* The 
total time for these operations is thus 1 h, 0 3 min. 


*New software from the plotter manufacturer should allow us to reduce this 
time . 
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C . Test Program 

It is very useful to have information on the quality of the spectrum as soon as 
possible after the recording of the inte rfe rogram. For this it suffices to cal- 
culate, at maximum resolution, three narrow regions (8K points wide, for 
instance) in three different spectral intervals; one in a region of maximum 
signal, within which are narrow lines, and the two others at the wings of the 
spectrum, where the signal is zero; one has, therefore, a simple means of 
determining the effective noise level and the signal/noise ratio in the spec- 
trum. The calculation of the FT is done from a secondary inte rfe rogram 
sampled at 8K equidistant points spaced at 1024K/8 intervals generated by a 
suitable convolution with the 10 primary samples (Ref. 8). Table 2 gives 
the corresponding calculation times. 

Thus, using a test taking 3 X 96 s =4 min 48 s, one has enough infor- 
mation about the spectrum to decide whether to perform the complete 
calculation . * 


*One can, equally well, calculate a complete spectrum of 10Z4K in 16 slices 
of 64K if one has access to an "Array Processor, " a specialized device for 
performing convolutions which works in parallel to the Central Processing 
Unit (CPU). The time for calculating the secondary inte rfe rogram then 
becomes 50 s, whence one may reduce an entire slice of 64K points in 50 + 
33.84 - 89 s and the entire spectrum in Z4 min. The Array Processor also 
reduces the time to interpolate the secondary spectral points from the 106 
primary points from 22 min to 40 s. These trials were performed on an 
Array Processor connected to the 360/ 65 at the French National Institute of 
Astronomy and Geophysics. 
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Table 1. Calculation times for various machines and programs. The 
number of samples in the interferogram from a path differ- 
ence of 0 to A is N. 

max 

NOTE ; ORDINATEUR = computer 

PLACE EN MEM. CENTR. = core used 


Ordinateur 


Programmes 


Fonction 


Place en 
m6m. centr. 


N= 0,5 K 


N — 1 K 


N = 


N — 

8 K 

N = 

16K 

N — 

32 K 

N = 

64 K 

N = 

128 K 

N = 

256 K 

iV = 

512 K 


classique 


rtelle 


4 N 


classique 


rdelle 


4 AT 


FFT 4 K \ avec bandes 


complexe 


5 N 


complexe 


SN 


360/75 

360/75 

360/75 

256* 

256 K 

256* 

FFT 32 

FFT 64 AT 

TF 1 000 A" 

N< 32 K 

N< 64* 

avec disques 

complexe 

r6elle 

rfelle 

5 N 

SN/2 

SN/32 

i 1 


3 mn 30 s 


14 mn 


55 mn 55 s 


3 h 43 mn 


14 h 56 mn 


78 s 


5 mn 14 s 15,33 s 


20 mn 58 s 32,45 s 


1 h 23 mn 68,41 s 


5 h 35 mn 


5 

mn 

15s 

9 

mn 

30 s 

21 

mn 


50 

mn 


1 h 45 

mn 


3 h 32 

mn 



0,30 s 


0,66 s 


1,34 s 


2,88 s 


5,85 s 


12,22 s 


25,45 s 


0,22 s 


0,44 s 


0,90 s 


1,85 s 


3,82 s 


7,86 s 


16,00 s 


33,44 s 


42,92 

1 mn 

24 s 

2 mn 

46 s 

5 mn 

27 s 

11 mn 

22 mn 

08 s 
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Table 2 


Calculation times for slices of 8 K, 32 K, and 64 K spectrum- 
primary points starting from 10 interferogram samples. 

NOTE; TF = FT = Fourier Transform 



8 K 

32 K 

64 K 

Calcul de rinterferogramme se- 
condaire 

82 s 

166 s 

265 s 

Calcul de la TF 

3,82 s 

16 s 

33,44 s 

Interpolation dans le Spectre. . . 

10 s 

40s 

80 s 

Temps total 

95,82 s 

222 s 

378,44 s 
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Fig. 13. The actual spectrum occuping the interval Act between cr ^ and o^. 

In the calculation from the samples 1(A), not only is the spectrum 
found, but also its images (aliases). 



Fig. 14. Plot of the calculation time for various machines and programs. 

The crosses are measured values. Dotted lines indicate extrapo- 
lations. 

Authors erratum ; In order to agree with the text and Table 1 
(which contains the most recent results), the curve for the 1000 K 
transform on a 360/75 should be lowered by a factor of 2. 
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V. PRELIMINARY RESULTS 


A . De scription 

Three almost identical interferometers are under construction: one for 
astronomical spectroscopy, one for laboratory absorption spectroscopy, and 
one for emission spectra in the laboratory. Only this latter has produced 
results, to date. The properties of very high resolution and precise fre- 
quency measures for this class of instrumentation having been demonstrated 
(Ref. 11), we now wish to explore the possibility of detecting very weak lines 
over a wide spectral range and in the presence of a considerable total energy. 

The sources used for the study of rare-earth spectra are magnetron- 
excited electrodeless discharge tubes; they have been difficult sources for 
Fourier spectroscopy until now. They are frequently very unstable and are 
particularly prone to go into a low-frequency oscillation at frequencies simi- 
lar to the sampling-rate and the internal modulation. Under such conditions, 
the fundamental noise limit of the detectors themselves has yet to be reached. 
Appreciable improvement can be made by a constant -intensity servo-control 
of the source itself. 

Figures 15, 16, 17, and 18 show various portions of a spectrum of 
thorium generated from the same interf erogram . The recording time used 
was 10 h (the minimum possible would have been 6 h) the step-length between 
samples was 0.455 pm (corresponding to p = 13). Consequently, the free 
spectral range Ac r goes from 0.91 pm to infinity (0 to 1 1,000 cm ^). The 
occupied spectral range is limited by the transparency of the silica envelope 
of the source (and the response of the PbS detectors) to 3.3 pm (3000 cm ^). 

A filter eliminates wavelengths less than 0.8 pm. The lines contained 
in the interval 0.8 to 0.91 pm (12,500 to 11,000 cm ^ ) appear in the output 
in the region 1 1,000 to 9,500 cm ; but as a consequence of the properties 
of a sine transform they appear negative and are very easy to distinguish 
from the others on high-resolution plots. An absorption spectrum, by con- 
trast, could not be recovered in this manner. 

The peak-to-peak amplitude of the internal modulation is equal to 
0 .70 pm (q = 20) . The efficiency is unity at \ = 1 ,70 pm (6,150cm ^ ) and is 
greater than 0 . 5 from 2, 400 cm ^ (4.2 pm) to 12, 000 cm ^ (0.83 pm) . 
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To summarize, the spectral region 0.8 to 3.3 fim is covered in a single 
recording at an acceptable efficiency and with no ambiguity in the interpre- 
tation. Absolutely no prior knowledge of the spectrum is needed. The PbS 
detector is not the best possible for use below 1.1 |um, but a useful output in 
a region already studied (in part by spectrographs) is thereby obtained. 

4 

Figure 15 shows the entire spectral region computed from the first 10 
points of the interferogram; thus the resolution is a factor of 100 worse than 
the final value. It allows an appreciation of the density of the spectrum. It 
should also be noticed that the continuous background has atmospheric water 
vapor absorption bands* near 3700 cm and 5400 cm superposed on it. 

The background reaches a maximum at 3660 cm A fraction of the emis- 
sion lines are due to iodine liberated by decomposition of the thorium iodide 
used in the source tube. 

- 1 5 

Figure 16 shows a slice of 400 cm computed from the first 10 points 
of the interferogram, and finally 40 cm ^ computed using all the interfero- 
gram (10 6 points). The maximum path difference attained was 45 cm*, with 

- 3 - 1 

the particular apodizing function used, the re solution was 6c r = 20 X 10 cm 

-3 - 1 

The line-widths measured in the spectra varied from 40 X 10 cm near 

-1 -3 -1 -1 

10,000 cm to 25 X 10 cm near 4,000 cm 

Figure 17 shows the profile of an isolated line near 4200 cm ^ , and 
Fig. 18, a very weak line near 3439 cm ^ (2.93 fim) superposed on a back- 
ground which shows H^O absorption (and noise). 

The reproducibility of weak lines is demonstrated in Fig. 16 by com- 
parison of two separate recordings. The differences in intensity are due to 
differences in the operating mode of the source, which is not precisely 
reproducible. All the plots are automatically normalized with respect to 
the strongest line in the spectrum (8, 904 cm ^). With respect to this line, 

the measured effective noise level (in a region of zero signal) B - - = 0 . 8 X 
«« 4 c II 

10 ; from the measure of the total energy in the spectrum, it is deduced 

that Q = 10 7 . 


*It is planned, eventually, to operate the system in a vacuum. [Translator's 
note: this was completed by July, 1971. J 
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Figure 19 shows results obtained with holmium; the recording 

conditions (N, p, q, Act, 6cr)were identical. A piece of the spectrum at 

4 

low resolution (10 points) is again pre sented; the high- re solution spectrum 

shows the characteristic eight-component structure for holmium (spin 7/2). 

7 

The measured value of the Q-factor is 2.5 X 10 ; the highest value obtained 
previously was 3 X 10^ (Ref. 2). 

B . Comparison With Results Obtained by Classical Techniques 

The study of rich line spectra is relatively easy in the visible region 
where a large grating spectrograph (Paschen-Runge or crossed-disper sion 
echelle grating) allows the entire spectrum to be recorded in a single expo- 
sure. Many lines of thorium, on the other hand, have been measured 
photographically with Fabry-Perot etalons and make very useful secondary 
wavelength standards (precisions of the order of 10 cm ) (Ref. 16). 

In the infrared, the only technique extensively used to date has been by 
scanning with a grating spectrometer. Steers (Ref. 17) recorded the thorium 

spectrum from 1 to 2.5 pm with an echelle grating; about 1000 lines were 

-2 - 1 

determined; the measured wave number errors varied from ±2 X 10 cm 
to ±1 X 10 ^ cm * 

A scanning Fabry-Perot interferometer gives much better precision; 
it has been used for noble gas spectra (Ref. 18). Since a separate recording 
is needed for each line, the method is slow and has not yet been employed 
for rare earths . 

At the Aime Cotton Laboratory, rare earth spectra are systematically 
explored with two grating interference spectrometers of the SISAM variety 
(Ref. 15), with which the speed gain over classical grating spectrometers 
is already significant. The technique and some results have been discussed 
by J. Verges (Refs. 19, 20, 21, and 22); the accent has been more on high 
resolution for the analysis of hyperfine structure, isotope shifts, and 
Zeeman effect rather than precise absolute positions. The method goes in 
two successive stages: 

(1) A preliminary scan of the entire spectrum is performed at 

0.2 cm ^ resolution using a first SISAM; the region covered is 
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4,000 to 12,000 cm< ^ . The effective recording time is of the 
order of 150 h and the operation takes about 3 weeks.* 

(2) After a first look at the spectrum, in the course of which the 

interesting lines are noted, an analysis at much higher resolu- 
tion (0.03 cm 1) is performed in very narrow spectral regions. 
The effective recording time for one structure is of the order of 
10 to 30 min; as a consequence of necessary adjustments and 
calibrations, the number of structures analyzed does not exceed 
10 to 15 per day. The study is necessarily limited to the strong- 
est lines only (50 in the case of holmium, for example). Fig- 
ure 20 compares the weakest line satisfactorily recorded with 
the SISAM to the result obtained by Fourier Transform spectro- 
scopy. The intensity of this structure is about 25 times greater 
than the weakest structure discernible in Fig. 19 (at 8152 cm ^ ) . 


*A similar study at longer wavelengths is performed by C. Morillon using 
a Girard grille spectrometer (Ref. 23). 
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INT ENSUE 1 



NOTE ; INTENSITE 1 = unit intensity, the peak to which this and all subsequent traces are 

normalized. 




Fig. 16. Spectrum of thorium (same interfer ogram as for Fig. 15) 

On Top: the spectral region indicated in Fig. 15, calculated from 

105 points (8cr = 0. 2 cm™^). 

Middle: 40 cm“^ slice calculated from 10^ points (5cr = 0.02 cm“^) 

Bottom: Same region and abscissas, but the vertical scale multi- 

plied by 10. A second spectrum recorded under the same 
conditions is superposed on the first. On this scale, the 
entire spectrum would be 80 meters long. 

NOTE; TRACES INFER IEURES = lower traces 
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Fig. 17. Spectrum of thorium 

Above: a region indicated on Fig. 15, calculated from 10^ points 

(6cr = 0. 2 cm”-*-) 

Below: a single line calculated from 10^ points (Scr = 0.0Z cm - ^) showing 

essentially a Lorentz profile. 
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THORIUM 


Below: 


Above: 



Fig. 18. Spectrum of thorium 

5 

the region indicated in Fig. 15, calculated from 10 points 
(8 or = 0.2 cm 

a very weak line calculated from 10^ points and showing the smallest 
width observed (0.025 cm ^). 
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4006 4007 cm 


Fig. 20. Spectrum of holmium 

Above: the weakest hyperfine structure of holmium recorded with a SISAM 
(6 cr = 0.03 cm“l, recording time 30 minutes) 

Below: same structure recorded by Fourier spectroscopy (same interfero- 

gram as for Fig. 19; 6c r = 0.02 cm ^). The entire spectrum is 320 
m long. 
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VI. CONCLUSION 


These first results have established the possibility of attacking 
problems needing high re solution and very wide spectral range using Fourier 
spectroscopy. In the case of complex emission spectra, the information 
contained in one 10^-point inte rfe rogram surpasses both in quantity and qual- 
ity that produced by existing methods in the course of research programs 
lasting several months. The calculated spectrum is available on a magne- 
tic tape in a format directly usable by the computer, and programs exist 
that automatically detect and measure the position and intensity of lines 
(Ref. 8). 

The next stage should be to undertake the measurement of absolute 
frequencies, which will require that the interferometer control system oper- 
ate from an international standard Krypton lamp. The interferometer will, 
eventually, be tried in the visible; its performance must then be compared 
to that of a large crossed-grating spectrograph with a Fabry-Perot etalon, 
compared to which there is no question of being able to obtain a significant 
speed gain. An increase in the number of samples to N = 5 X 10^ seems 
conceivable without modification of the instrument; spectroscopic problems 
justifying such an expansion are doubtless rare. 

A real-time digital compute r, specially adapted to the interferometer, 
is also in the process of construction.* It uses the classical algorithm (that 
of Cooley and Tukey not being usable in a real-time calculation). The basic 
cycle -time is 1 pm s/input point/output point; thus it is, for this particular 
problem, 1000 times faster than the 7040 (Table 1, column 2, and Fig. 14, 
the curve labeled n 7040 classique") and about 100 times faster than the 
360/75. It can calculate up to 20, 000 points in a spectrum starting from an 
inte rfe rogram recorded at a rate of 50 points/ second, without any limitation 
on N . In itself sufficient for all problems in which M < 20, 000, in other 
cases it allows a portion of the spectrum to be examined in real time. 


^Translator 1 s Note: This device was successfully completed early in 1971. 
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